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Abstract 

The proton-to-pion ratios measured in the BRAHMS experiment for Au+Au and p+p collisions at -yTvJviv = 62.4 and 
200 GeV are presented as a function of transverse momentum and collision centrality at selected pseudorapidities in 
the range of to 3.8. A strong pseudorapidity dependence of these ratios is observed. We also compare the magnitude 
and /^-dependence of the p/n + ratios measured in Au+Au collisions at yfsNN = 200 GeV and rj * 2.2 with the same 
ratio measured at y/s^^ = 62.4 GeV and rj = 0. The great similarity found between these ratios throughout the whole 
Pt range (up to 2.2 GeV/c) is consistent with particle ratios in A+A collisions being described with grand-canonical 
distributions characterized by the baryo-chemical potential pn. At the collision energy of 62.4 GeV, we have observed 
a unique point in pseudorapidity, rj = 3.2, where the p/n + ratio is independent of the collision system size in a wide 
/77-range of 0.3 < pr < 1.8 GeV/c. 
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1. Introduction 

The ongoing dialogue between theory and experi- 
ment is delivering an increasingly clear picture of the 
QCD phase diagram in its partonic and hadronic phases. 
In particular, recent mid-rapidity RHIC results are con- 
sidered as evidence for a smooth cross over from a dense 
and opaque partonic phase into a hadronic gas at tem- 
peratures around 170 MeV and small values of baryo- 
chemical potential pg 01. Measurements of the rela- 
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tive abundances of hadronic species constrain statisti- 
cal models used to describe the chemical freeze-out in 
nucleus-nucleus interactions at different energies. Some 
of these models show a remarkable behavior at low 
baryo-chemical potential where the curves of temper- 
ature vs. baryo-chemical potential representing chem- 
ical freeze-out tend to merge with the phase boundary 
between partonic and hadronic media yfl. In such a pic- 
ture, hadrons are produced very close to freeze-out, and 
one could safely state that the features of the partonic 
medium are transmitted to the final bulk hadrons via 
hadronization processes. Indeed, such appears to be the 
case for the constituent quark scaling found in elliptic 
flow V2 measurements J3] as well as the enhancement 
of baryon-to-meson ratios that scale with the collision 
size around mid-rapidity (low /is) This increase 

in the baryon-to-meson ratios (« 1 for the so-called in- 
termediate pr values ranging from 2 to 6 GeV/c), first 
reported by the PHENIX Collaboration fl, deviates 
remarkably from calculations that include parton frag- 
mentation in the vacuum. The /^-dependence of the 
baryon-to-meson ratio in the intermediate pj range ap- 
pears to be related to a modification of hadronization in 
the partonic medium created by the collision. This ef- 
fect could be caused by the different quark content of 
the baryons and mesons Q] or because of their different 
masses through the effect of radial flow @|. 

Both radial flow and in-medium quark coalescence 
are expected to enhance protons over pions at interme- 
diate pj. In particular, the PHENIX pjif data at mid- 
rapidity is well described by the Greco, Ko, and Levai 
quark coalescence model |9J] where the introduced co- 
alescence involves partons from the medium (thermal) 
and partons from mini-jets. The Hwa and Yang quark 
recombination model I Kill describes well the BRAHMS 
and PHENIX p/n + ratios around mid-rapidity. These 
mid-rapidity results are consistent with hadronization in 
the intermediate pj range being dominated by parton re- 
combination with negligible final state interactions be- 
tween the produced hadrons. In contrast, at forward ra- 
pidity (large a significant gap between the temper- 
ature of the transition from the partonic to the hadronic 
phase, T c , and the temperature of chemical freeze-out is 
predicted by QCD lattice calculation fiUl . In that envi- 
ronment, hadronic re-scattering may play an important 
role, making statistical models more suitable to describe 
particle abundances at high rapidity 1112111311 . 

We present in this Letter the p/n ratios for both 
charges extracted from Au+Au collisions at different 
centralities as well as p+p at y/s^N = 62.4 and 200 GeV. 
The ratios are presented as function of pseudorapidity 
and transverse momentum pj. The system size depen- 



dence of the ratios is investigated and compared to the 
smallest system; p+p collisions. The comparison of 
the p/n + measured at two different energies and pseu- 
dorapidities is used to establish a connection between 
particle abundances in these systems and their possible 
description by a grand-canonical distribution. We also 
present an interesting feature in the pseudorapidity de- 
pendence of the p/n + ratio at -\/snn = 62.4 GeV. There 
the ratio has the same magnitude and pj dependence 
at T] = 3.2 when measured at all Au+Au centralities as 
well as p+p collisions. We compare the ratios extracted 
from the most central Au+Au events to a hydrodynami- 
cal model and to a partonic cascade model that includes 
hadronic re-scattering in the final state. 

2. Experimental setup 



The BRAHMS detector setup yj] consists of two 
movable, small acceptance spectrometers: the Mid- 
Rapidity Spectrometer (MRS), which operates in the 
polar angle interval from 90° > 6 > 30° (correspond- 
ing to a pseudorapidity interval of < 77 < 1.3) and 
the Forward Spectrometer (FS), which operates in the 
range from 15° > 6 > 23° (2 < 77 < 4). The 
BRAHMS setup also includes detectors used to deter- 
mine global features of the collision such as the over- 
all charged particle multiplicity, the collision vertex and 
the centrality of the collision. The MRS is composed 
of a single dipole magnet placed between two Time 
Projection Chambers (TPC), which provide a momen- 
tum measurement. Particle identification (PID) is based 
on Time-of-Flight (TOF) measurements H 1 411 - The FS 
has two TPCs, which are capable of track recognition 
in a high multiplicity environment close to the inter- 
action region, and, at the far end of the spectrometer, 
three Drift Chambers. In the aggregate, the FS can de- 
tect particle track segments with high momentum res- 
olution (Sp/p = 0.0008/? at the highest field setting) 
using three dipole magnets. Particle identification in the 
FS is provided by TOF measurements in two separate 
hodoscopes for low and medium particle momenta, re- 
spectively. High momentum particles are identified us- 
ing a Ring Imaging Cherenkov detector (RICH) fl5ll . 

3. The analysis 

The analysis reported in this Letter consists of the 
comparison of proton and pion yields as a function of 
pj for several pseudorapidity intervals. The analysis is 
carried out in 77 versus pj space, since for any given an- 
gle and field configuration the acceptance is the same 



for pions and protons. For a given r\-pj bin the p/rc 
ratios are calculated on a setting by setting basis. In or- 
der to avoid mixing different PID techniques, which can 
lead to different systematic uncertainties, the ratios are 
calculated separately for PID using TOF and the RICH 
detectors. In this way, all factors such as acceptance 
corrections, tracking efficiencies, trigger normalization 
and bias related to the centrality cut cancel out in the 
ratio. The remaining species dependent corrections are: 

(i) decays in flight and interactions with the beam 
pipe and other material, and 

(ii) the PID efficiency correction. 

The corrections for (H} are determined from the sin- 
gle particle response of pions and protons in a realistic 
GEANT H3] model description of the BRAHMS ex- 
perimental setup. The magnitude of this correction on 
the particle ratios depends on the particle momenta and 
the spectrometer positions, but does not exceed 6%. We 
estimate that the overall systematic uncertainty related 
to this correction is at the level of 2%. 

The TOF PID is done separately for small momen- 
tum bins by fitting a multi-Gaussian function to the ex- 
perimental squared mass m 2 distribution and applying a 
±3cr cut to select a given particle type. For measure- 
ments done with the FS spectrometer in the momentum 
range where pions overlap with kaons (usually above 
3.5 GeV/c), the RICH detector can be used in veto mode 
to select kaons with momentum smaller than the kaon 
Cherenkov threshold, which is about 9 GeV/c. This pro- 
cedure leads to a relatively clean sample of pions with 
some contamination by kaons having spurious rings as- 
sociated in the RICH counter. Together with the kaon 
- proton overlap at larger momenta, this contamination 
effect is a source of systematic errors which have been 
estimated to be in the order of a few percent for all pj 
values, except for low pjin the pjii~ ratios, where these 
uncertainties can reach high values (« 15%). At mid- 
rapidity the systematic uncertainty reaches a value of 
5% at pj > 2.5 GeV/c due to the limited kaon to pion 
separation at large momenta. 

The RICH PID is also based on the particle separation 
in the rrr versus momentum space. The RICH provides 
direct proton identification above the proton threshold 
momentum which is about 15 GeV/c. However, an ad- 
ditional proton identification scheme is possible below 
this value but above the kaon threshold. In this momen- 
tum range a proton is associated with tracks having mo- 
menta above the kaon threshold, but no RICH signal 
(veto mode). The veto proton yields are corrected for 
pion and kaon contamination due to a small RICH in- 
efficiency. The RICH inefficiency was determined by 



studying yields of tracks identified in the TOF detec- 
tor as pions but having no associated ring in the RICH. 
This study was done for a low momentum range with a 
good kaon/pion separation in TOF. It is found that the 
pion efficiency grows rapidly from the pion threshold 
(~ 2.3 GeV/c) and reaches a constant value of about 
97% above 4 GeV/c. The RICH inefficiency found for 
pions is then applied to kaons (and the other species) 
assuming that Cherenkov radiation depends only on the 
y factor of the particle. A more detailed description 
of the RICH inefficiency analysis is given in IU5I1 and 
a forthcoming publication. There are two sources of 
systematic uncertainties associated with the RICH PID, 
namely, the uncertainty on the RICH inefficiency, esti- 
mated to be at the level of 10%, and the overlap in m 2 
between pions and kaons having momenta above about 
30 GeV/c. 

No corrections were applied for weak decays. How- 
ever we apply cuts on track and event vertex matching 
to limit the effects of feed down on particle yields. The 
ranges of the vertex cuts are determined by the mea- 
sured spatial resolution of the particle track projection to 
the event vertex, such that 97% of primary tracks are in- 
cluded. Using AMPT lfl7ll model calculations we found 
that these vertex cuts lead to remnant contamination of 
proton and pion yields mainly due to A hyperons and 
K° meson decays. At mid-rapidity we estimated that 
this contamination leads to a 10%(7%) enhancement 
of p/n + (p/n~) at pj = 3.0 GeV/c. This enhancement 
increases toward low transverse momenta and reaches 
20%(14%) for p/tt + (p/tt-) at p T = 1.5 GeV/c. The 
level of contamination decreases gradually towards for- 
ward rapidities and at rj « 3 the enhancement reaches 
values of 6%(5%) for pln+ipliC) at p T = 3.0 GeV/c 
and 10%(7%) for p/n + (p/n-) at p T = 1.5 GeV/c. We 
also found from the model calculations that the ex- 
tracted contaminations are proportional to the primary 
A/p ratios (A/p ratios for negative species). Thus if the 
AMPT model under predicts the data by 20% then the 
respective corrections due to contaminations should be 
increased by approximately 20%. 

4. Results 

Figure [TJ shows pln + ratios measured in Au+Au col- 
lisions at -\/snn - 200 GeV in two centrality classes 
of events, namely, - 10% (solid dots) and 40 - 80% 
(open squares). The centrality selection is based on 
the charged-particle multiplicity measured in the range 
-2.2 < 77 < 2.2 as described in [18]. The shaded boxes, 
plotted for the most central events, represent the sys- 
tematic uncertainties discussed in the previous section. 
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Figure 1: Centrality dependent p/n + ratios for Au+Au system collid- 
ing at ys,vjv = 200 GeV for central (0 - 10%) and semi-peripheral 
(40 - 80%) reactions in comparison with p+p collision data at the 
same energy. The vertical bars represent the statistical errors and the 
shaded boxes (plotted only for central Au+Au) show the systematic 
uncertainties. 



However, they do not include the uncertainties related 
to the weak decay contamination. The ratios extracted 
from p+p data at the same energy are plotted for com- 
parison (solid triangles). The pj coverage depends on 
the pseudorapidity bins and for central Au+Au colli- 
sions extends up to pj — 4 GeV/c for rj = 2.6. At low 
Pt (< 1.5 GeV/c) the p/n + ratios exhibit a rising trend 
with a weak dependence on centrality. A significant de- 
pendence on centrality begins above 1.5 GeV/c. The 
ratios appear to reach a maximum value at pj around 
2.5 GeV/c wherever the pj coverage is sufficient to de- 
termine this limit. The maxima of the ratios are greater 
for more central events and, at i] = 3.1, are equal to 
about 2.5 and 1.5 for the - 10% and 40 - 80% cen- 
trality bins, respectively. The p+p ratios are consistent 
with Au+Au data at low pj and begin to deviate sig- 
nificantly above pr = 1 GeV/c. At 77 = 3.1 a maxi- 
mum value of the ratio of 0.55 is reached in p+p col- 
lisions which is a factor of 4.5 smaller than that ob- 
served for central Au+Au reactions. We have performed 
PYTHIA II 1911 calculations for nucleon-nucleon interac- 
tions in each possible isospin state. The calculations of 



Figure 2: The ratios of p/7r~measured in Au+Au and p+p collisions 
at y.SAw = 200 GeV, are plotted for the same centrality samples used 
in Fig.[T] 



the p/tt ratio vs. pj show no significant dependence on 
the isospin of the initial NN system. Thus the difference 
between p+p and Au+Au is not an isospin effect. 

The values of the p/n~ ratios plotted in Fig. [2] are 
significantly lower than the p/n + ratios (note the differ- 
ence in the vertical scale), however, the centrality de- 
pendence shows the same features as those observed 
in the p/n + ratios, namely, that the ratios for different 
centrality classes are consistent with each other up to 
Pt ~ 1.2 GeV/c and a strong dependence on centrality 
appears at larger transverse momenta reaching a max- 
imum at similar pj as the positive particles. Looking 
at the p+p data alone, one notes the difference in shape 
between the p/n + and pin' ratios: a clear shift of the 
pin' peaks towards lower pj, as well as a much broader 
p/n + peaks. The large difference between the Au+Au 
and p+p both in shape and overall magnitude may re- 
flect significant medium effects in Au+Au collisions at 
yfsim = 200 GeV in the pseudorapidity intervals cov- 
ered. 

In Fig. [3]we explore the possible connection between 
the measured p/n ratios in extended systems and the 
baryo-chemical potential pg used to characterize them 
as statistical systems. Such a connection is done by 
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Figure 3: p/n + ratio in Au+Au collisions for rj = 0.0 at ys/vw = 
62.4 GeV marked with open triangles and the Au+Au reactions for 
r\ = 2.2 at ysjviv = 200 GeV marked with the black triangles. Stars 
shown reference mid-rapidity data for p+p at ys/viv = 62.4 GeV. 



comparing two rapidity ranges that have the same pip 
ratio but different yfsim- These are r\ — at 62 GeV 
and rj = 2.2 at 200 GeV. The pseudorapidity intervals 
selected for this comparison correspond to similar ob- 
served pip ratios of approximately 0.45. The very sim- 
ilar pip ratios in these two systems has been attributed 
to their having a common value of the baryo-chemical 
potential p B of ~ 62 MeV BED]. The similarity of 
proton-to-pion ratios for these selected heavy ion colli- 
sions suggests that the baryon and meson production at 
the pr interval studied (up to 2 GeV/c) is dominated by 
medium effects and is determined by the bulk medium 
properties. The considerably lower values of the p/n + 
ratio measured in the p+p system at -\fs = 200 GeV, 
shown with stars in Fig. [3] can also be construed as 
strong indication that medium effects are the source of 
the observed enhancement of the p/n as function of pj 
in the nucleus-nucleus collisions. 

The three panels of Fig. [4] display the p/n + ratio ex- 
tracted from p+p and Au+Au collisions at yijviv = 
62.4 GeV measured at 3 different pseudorapidity bins; 
77 = 2.67 (top panel), rj = 3.2 (middle panel) and rj 
= 3.5 (bottom panel). At 77 = 2.67, p/ti + is greater 
in Au+Au than in p+p collisions while at rj — 3.5 
the situation is reversed. The middle panel of Fig. [4] 
shows the p/n + ratios from Au+Au and p+p collisions 
at ysiviv = 62.4 GeV at 77 = 3.2. This "crossing point" in 
pseudorapidity shows a remarkably complete overlap of 
the ratios as function of pj not only for p+p and 0- 10% 
central Au+Au reactions, but also for other Au+Au cen- 
tralities, namely 10-20%, 20-40% and 40-80%. Such 
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Figure 4: p/rr + ratio from p+p and Au+Au collisions at yiww = 
62.4 GeV for rj = 2.67 (top panel), r\ = 3.2 (middle panel) and r\ = 
3.5 (bottom panel). The curve drawn in the middle panel is a polyno- 
mial fit to the most central events, its purpose is to guide the eye and 
bridge the gap in pj . 



a universal shape for the p/n + ratios implies, that for 
each centrality at this pseudorapidity, the nuclear modi- 
fication factors for protons and pions are consistent with 
each other at all measured values of pj. The observed 
crossing of p/n + ratios for different sizes of the collid- 
ing systems is consistent with recent BRAHMS results 
on the centrality dependence of net-proton rapidity dis- 
tribution in Au+Au reactions at yfsim = 200 GeV B22I1 . 
The data show an increased baryon transport towards 
mid-rapidity in central collisions. Such increased stop- 
ping power dissipates the beam energy to form a denser 
system where recombination mechanisms favor proton 
production at intermediate values of pj. Both effects 
will produce higher values of the p/n + ratio in larger 
systems. On the other hand, at very forward rapidities 
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Figure 5: Rapidity evolution of p/n + and p/7r for - 10% central Au+Au reaction at -yjsp/N = 200 GeV (solid points), and calculations from the 
THERMINATOR (dashed line) and the AMPT models (solid line). 



(around one unit below the beam rapidity) a bigger frac- 
tion of the measured particles are protons and their num- 
bers are even higher for lighter systems due to reduced 
stopping power. This observation is consistent with the 
reversed p/n + dependence on the system size seen in 
the bottom panel. A crossing of p/n + yields ratio at 
the energyof 62.4 GeV is predicted by UrQMD El . 
HIJING IB and AMPT model calculations. However, 
these models predict the location of the crossing point 
in the interval 2 < y < 2.5, which is almost one unit 
of rapidity lower than the observed value. These exper- 
imental results then provide a strong constraint on the 
theoretical description of baryon number transport and 
associated energy dissipation in relativistic nuclear re- 
actions. It is worth noting that at 200 GeV (see Fig. [TJ, 
even for 77 = 3.8, the p/n + ratios measured in Au+Au re- 
actions are larger than those observed in p+p collisions, 
so the possible crossing point at high energy is located 
at larger rapidity, beyond the experimental acceptance. 

5. Model comparisons 

To interpret these results, theoretical models of 
nucleus-nucleus collisions are confronted with the data. 
The p/tt + ratio measured in central Au+Au collisions 
at -\fsNN = 200 GeV, shown on the top row of Fig. |5j 
has a strong growth with rapidity; starting from a value 
of about 1.0 at 77 « and pj « 3 GeV/c and reach- 



ing a value of 2.5 at 77 ~ 3.5 and pj « 3 GeV/c. In 
contrast, the p/n~ ratio (bottom row) decreases with in- 
creasing rapidity ( from * 1 at 77 * to 0.4 at 77 =s 3.5). 
Note the different vertical scales for positive and nega- 
tive charged particles. Figure[5]compares our results for 
- 10% central Au+Au reactions at y/s NN = 200 GeV 
to two model calculations based on: the THERMINA- 



TOR model |25J1 (black dashed lines) and the AMPT 
model (solid lines). THERMINATOR is a 1 + 1-D hy- 
drodynamic model that incorporates statistical particle 
production (including excited states), which in turn, has 
chemical potentials with parameterized rapidity depen- 
dence. Alternatively, the AMPT (A Multi-Phase Parton 
Transport) model includes mini-jet parton production, 
parton dynamics, hadronization according to the LUND 
string fragmentation model, and final state hadronic in- 
teractions in determining the final particle production. 
In these simulations we duplicate experimental condi- 
tions regarding the particle contamination due to weak 
decays, the tracks were required to point to the event 
vertex under the same conditions as those applied to the 
experimental data. 

The THERMINATOR model tracks well the trends 
of the data up to pj ~ 2.5 GeV/c for both positive 
and negative particles. At 77 = this model describes 
well the experimental p/ji + and pjii~ ratios in the in- 
termediate pj range covered by the data (1 < pj < 3 
GeV/c), however, it predicts the p/n~ ratios with val- 



ues well above those measured at forward rapidities. 
At large p T (> 3 GeV/c) THERMINATOR fails to de- 
scribe the data. The mismatch is clearly seen in the 
rj = 3.05 plot of Fig. [5] where the data have the best 
pj coverage. This mismatch is attributed to the fact 
that the model does not include the production of jets. 
The AMPT model can qualitatively describe the pseu- 
dorapidity trends, but fails to quantitatively describe the 
data, namely, the model under predicts p/n + and over 
predicts pjif ratios. We note that the AMPT describes 
the p/n + ratios reasonably well for semi-peripheral re- 
actions (not shown). 

6. Summary 

We present the /^-dependence of the p/n ratios mea- 
sured in Au+Au and p+p collisions at energies of 62.4 
and 200 GeV as a function of pseudorapidity and col- 
lision centrality. The data provide the opportunity to 
study baryon-to-meson production over a wide range of 
the baryo-chemical potential /ub- For Au+Au collisions 
at y/sNN = 200 GeV the p/n + and pjif~ ratios show a 
noticeable dependence on centrality at intermediate pj 
with a rising trend from p+p to central Au+Au colli- 
sions. We find that p/n + ratios are remarkably similar 
for central Au+Au at r\ ss 2.2 at ysjv/v = 200 GeV, 
and central Au+Au at r\ * at ysMv = 62.4 GeV, 
where the bulk medium is characterized by the same 
value of pip. This observation, together with the ob- 
served centrality dependence suggests that particle pro- 
duction for intermediate pj values is governed by the 
size and the chemical properties of the created medium 
for the systems and pseudorapidity range studied. It is 
also shown that the statistical model describes well the 
mid-rapidity p/n + and p/n~ ratios for central Au+Au at 
yfsmi = 200 GeV, and it also tracks well the data trends 
at high rapidity. Finally, the Au+Au and p+p measure- 
ments at -\/snn = 62.4 GeV show that the p/n + ratios 
for p+p and for all analyzed Au+Au centralities cross 
simultaneously at the same 77 value (« 3.2) and are con- 
sistent with each other in the covered pj range, e.g., 
from 0.3 GeV/c up to 1.8 GeV/c. 
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